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The nature of the mid-infrared population from optical identifications of the 
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F. La Franca\ C. Gruppioni^, I. Matute\ F. Pozzi^-^, C. Lari'^, M. Mignoli^, G. Zamorani^, D.M. 
Alexander^, F. Cocchia^, L. Danese^, A. Franceschini^, P. Heraudeau^, J.K. Kotilainen^°, M.J.D. 
Linden- V0rnle^\ S. Oliver^^, M. Rowan-Robinson^^, S. Serjeant^"*, L. Spinoglio^^ & A. Verma^^ 

ABSTRACT 

We present a multi-wavelength catalog (15 /im, R, K-band, 1.4 GHz flux) plus spectroscopic 
identifications for 406 15 /im sources detected in the European Large Area ISO Survey (ELAIS) 
region SI, over the flux density range 0.5<S'i5^ni<150 niJy. 332 (^^82%) sources are optically 
identified down to R~23.0. Spectra or bona fide stellar identifications are obtained for 290 objects 
(~88% of the optically identified sources). The areal coverage, mid-infrared (MIR) and optical 
completeness of the sample are discussed in order to allow statistical and evolutionary analyses. 
Two main spectroscopic classes have been found to dominate the MIR extragalactic population: 
z < 0.5 star-forming galaxies (from absorbed to extreme starbursts: i^Ly(15/im)w 10^ — 10^^ L©), 
which account for ^75% of the sources, and Active Galactic Nuclei (AGN; both type 1 and 
2), which account for ^25% of the sources. About 20% of the extragalactic sources are dust- 
enshrouded starburst galaxies [e(a) spectra], and all the starburst galaxies appear more dust 
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extincted in the optical than nearby normal galaxies. We also identified 91 stellar objects (~22% 
of the MIR sources). The counts for star burst galaxies and AGN down to 0.6 mJy have been 
derived. A general trend is found in the optical-MIR spectral energy distribution (SED) of the 
galaxies, where the MIR-luminous objects have larger MIR to optical luminosity ratios. Based 
on a variety of analyses, wc suggest that the ELAIS sources fainter than R^23 arc luminous 
and ultra-luminous MIR galaxies (LIG-ULIGs; z/Lj^(15/zm)=10^^-10^^ Lq) at intermediate red- 
shifts (2;=0.5-1.5), and that consequently the present sample is virtually 100% spectroscopically 
complete up to z=0.5. 

Subject headings: Cosmology: observations; Galaxies: distances and redshift; Galaxies: evolu- 
tion; Galaxies: active; Galaxies: starburst; Infrared: galaxies 

1. Introduction 

The ISOCAM instrument (Cesarsky et al. 1996) on board of the Infrared Space Observatory (ISO; 
Kessler et al. 1996) has been able to unveil at 15 fim a population of MIR galaxies with two orders of 
magnitude fainter flux densities than observed by the Infrared Astronomical Satellite {IRAS}. These galaxies 
result to be ten times more numerous than expected if there were no evolution from 2: = up to z = 1 — 1.5 
(see Elbaz et al. 1999). These results are confirmed by the detection of a substantial diffuse cosmic Infrared 
Background in the 140 /im - 1 mm range (Puget et al. 1996; Hauser et al. 1998; Fixsen et al. 1998; Lagache 
et al. 1999), which suggests a strong evolution of the galaxies in the IR (stronger than observed at any 
other wavelengths). The suggested cvohition can either imply a larger fraction of galaxies experiencing an 
IR luminous phase in the past, or that MIR galaxies were more luminous in the past (see e.g. Chary & Elbaz 
2001). 

So far, the nature of the sources responsible for the strong evolution observed in the MIR band has 
only been studied in small fields, containing sources well known at other wavelengths: the Hubble Deep 
Field North (HDF-N; Ausscl ct al. 1999a,b) and South (HDF-S; Oliver et al. 2002; Mann et al. 2002; 
Franceschini et al. 2003), the Canada-France Redshift Survey (CFRS) 1452-1-52 field (Flores et al. 1999) 
and the ELAIS-S2 field (Pozzi et al. 2003). This has yielded a small but meaningful sample of sources 
sufficient for most multiwavelength studies. Most of these sources can easily be identified with relatively 
bright optical counterparts (/ < 22.5) with a median redshift of 2; ~ 0.7—0.8 imposed by the MIR fc-correction 
(Ausscl ct al. 1999a,b; Flores ct al. 1999). 

However, to obtain reliable source counts and statistically significant information about the nature 
and evolution of the MIR extragalactic source population, large-area surveys covering a wide range of flux 
densities are required. In particular, it is very important to bridge the gap between the IRAS surveys 
(S'ls^m > 200 mJy) and the Deep ISOCAM surveys (0.1 < Si^^m < 2 - 3 mJy). ELAIS (Oliver ct al. 
2000), with its broad flux range (0.5< S'15/zm <150 mJy) and several thousand sources detected over a ~12 
deg^ area (Rowan- Robinson et al. 2004), is the crucial survey to investigate the nature and evolution of the 
MIR sources in the region where the source counts start diverging from the non-evolution predictions (a few 
mJy at 15 /im; Elbaz et al. 1999; Gruppioni et al. 2002) and with data not affected by cosmic variance 
(which might be a strong effect in the narrow area surveys) . 

In order to better understand the nature of the MIR extragalactic population on the basis of a larger 
identification statistics than the ones obtained so far, we performed photometric and spectroscopic identifi- 
cations for a source sample selected at 15 fim in the ELAIS southern field SI (~4 deg^). Here we present the 
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results of the identification of 406 objects, which form a complete and highly reliable sub-sample drawn from 
the 15 /Ltm catalog of Lari et al. (2001). The areal coverage, MIR and optical completeness are discussed in 
detail to allow statistical and evolutionary analyses with this sample. The counts down to 0.6 mJy for each 
class of source have been derived. The multiwavelength (radio, MIR, near-infrared (NIR), optical) properties 
of the 15 /im sources arc investigated in order to study the nature and evolution of the infrared cxtragalactic 
population. The infrared luminosity function and the evolution of 15 /im galaxies and AGN are discussed 
in two companion papers (Pozzi et al. 2004; Matute et al. 2004, in preparation). A preliminary analysis of 
the AGN evolution has been presented in Matute et al. (2002) ^ 

This paper is structured as follows: in section 2 we describe the 15 /um sample; in section 3 we present the 

optical CCD observations, the source identification and the cross-correlation with catalogs in other bands; 
in section 4 we describe the spectroscopic observations, the classification and the spectroscopic catalog; in 
section 5 we present the counts; in section 6 we discuss the nature of the MIR galaxies. 

Throughout this paper we will assume Hq — 75 km s^^ Mpc^^, flm = 0.3 and f^A = 0.7. 

2. The ELAIS-Sl Mid-Infrared Sample 

The primary ISO sm-vcy in the ELAIS regions was restricted to two bands: 15 /um (using ISOGAM), 
and 90 /im (using ISOPHOT; Lcmke et al. 1996). Additional observations were performed in some parts of 
the survey area at 6.7 /im and 175 /^m. The main survey area, consisting of ~12 deg^, is divided into four 
fields distributed across the sky in order to decrease the biases due to cosmic variance: three in the northern 
hemisphere (Nl, N2 and N3) and one in the southern hemisphere (SI). Details of the survey can be found 
in Oliver et al. (2000). The final band merged catalog (6.7, 15, 90 and 175 /xm) of the whole ELAIS survey 
is in Rowan- Robinson et al. (2003)^. 

The southern field, SI, centered at a(2000) = OO'' 34" 44.4^ (5(2000) = -43° 28' 12", covers an area of 
the sky of about 2° x 2° . In this field, a new method of data reduction (LARI technique) , especially developed 
for the detection of faint sources, has been successfully applied, producing a catalog (complete at the 5a 
level) of 462 sources in the flux density range 0.5 150 m.Jy'^. Details about the data reduction technique 
and the source catalog have been presented in Lari et al. (2001)^, while the sample completeness and the 
source counts at 15 /um obtained from that catalog have been presented and discussed by Gruppioni et al. 
(2002). 

The whole SI area has been surveyed in the radio with the Australia Telescope Compact Array down to 
SiAGHz — 0.4 mJy (5(t, Gruppioni et al. 1999)^. The radio catalog in SI consists of 652 sources detected at 
the 5a level at 1.4 GHz. The radio- MIR correlation for starburst galaxies has been discussed by Gruppioni 
et al. (2003). About 40% the area has been observed by BeppoSAX down to a 2-10 keV sensitivity of ^10"^^ 
erg s~^ cm~^ by Alexander et al. (2001). 



'^A collection of papers and scientific results related to the ELAIS southern fields is available at 
http://www.fis.uniroma3.it/~ELAIS_S. 

^ See http: // astro. imperial. ac. uk / elais / . 

^After flux corrections as explained at the end of §2. 

"^The catalog is available at http://www.bo.astro.it/~elais/catalogs/ELAIS_CAM_15micron_Sl.TAB. 
^Available at http://www.bo.astro.it/~elais/catalogs/ELAIS_RADIO_Sl.TAB. 



-4- 



In this work we restricted our analysis to a highly reliable subsample of 406 sources out of a total 462 15 
/um sources. The ELAIS-Sl area has been observed by ISOCAM at 15 /xm with 9 separate rasters (see Figure 
1). 20 out of the 462 sources were excluded as they were detected twice at the edges of overlapping rasters 
(in this case the sources with the highest S/N ratio have been selected). Moreover, 36 sources were excluded 
as wc decided to use only those sources with the highest reliability (REL=0) after a visual inspection of 
their pixel histories. The ISOCAM observations of the central (Sl_5) raster have been repeated three times 
(instead of one) and thus in this part the survey reaches fainter fluxes. As a consequence, we have divided 
the sample into two regions: a) Sl_5 (i.e. raster 5) which covers 0.55 deg^ and reaches a 20% completeness 
at ^^0.7 mJy, and b) Sl_rest (i.e. all the remaining rasters) which covers 3.65 deg^ and reaches a 20% 
completeness at ~1.1 mJy. Note that in this work, as well as in the counts paper (Gruppioni et al. 2002), 
we used the 15 /xm fluxes of the ELAIS-Sl sources reported by Lari et al. (2001), corrected following the 
detailed recipe described by Gruppioni et al. (2002). The completeness functions of the two regions as a 
function of the corrected 15 fim flux arc shown in Figure 2 and are reported in Table 1, while the corrected 
15 /xm fluxes are reported in Table 2 and Table 3. 

3. Optical photometry and multi-wavelength identifications 

3.1. CCD Observations and Data Reduction 

In order to obtain R-band optical photometry of the ELAIS-Sl region down to magnitudes fainter 
than the typical Schmidt plates limits (e.g. COSMOS^), we have carried out a CCD imaging campaign 

at the European Southern Observatory (ESO^). The whole SI region has been covered by 168 overlapping 
CCD exposures obtained with DFOSC mounted on the 1.54 Danish/ESO telescope at La Silla (Chile). The 
observing runs were carried out during 12 nights (almost all photometric) during September 1996, October 
1996 and September 1997. The CCD had 2052x2052 pixel array, with a pixel scale of 0.39", and a field of 
view of 13.3x13.3 arc-minutes. The filter used was the Bessel R (ESO #452). About 97% the ELAIS-Sl 
area has been observed in the R-band (see Figure 1). 

Photometric calibration was carried out by frequently observing (every ~1-1.5 hours) the standard field 
stars "T Phe" (q(2000) = OO'' 27™ 49.0^ <5(2000) = -46° 48' 02"; Landolt 1992). Those fields (about 10%) 
observed without photometric conditions were later calibrated using shorter integration images taken during 
photometric weather. Flat fields were obtained from blank fields in the sky at sunset and sunrise every night. 
The integration time for all science frames was 20 minutes. Standard data reduction was carried out with 
the MIDAS and IRAF reduction packages. Every science frame was subtracted for bias, and corrected for 
flat field. 

Source searching on the final reduced images was carried out with SExtractor (Bertin & Arnouts 1996) 
The "MAG_BEST" source magnitudes have been used. The R-band magnitudes were computed after cor- 
rection for airmass extinction, and without a color term, i.e. assuming a V — R = color for all the objects. 
Due to the long observing campaign the characteristics of the CCD observations (e.g. seeing and magnitude 
limits) varied from field to field. The magnitude limits vary from 22.5 to 23.0 at 95% completeness level. 



®http: / /xip. nrl.navy.mil/www_rsearch/RS_form. html. 
''http://www.eso.org. 
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The astrometry has been computed by using the USN02 catalog* as reference, with an internal accuracy of 
0.26" (rms), and an absolute radial accuracy (compared to USN02) of 0.38" (rms). 

3.2. Identification of the 15 /xm Sources in the R band 

For the optical identifications of the ISOCAM sources we used the R magnitudes and positions given 
in the CCD catalog. Ton sources were not observed by the CCD imaging. For 7 out of these 10 sources 
the R-band USN02 magnitudes have been used. The remaining 3 sources (ELAISC15_J003559-430232, 
ELAISC15_J003912-442051, ELAISC15_J004023-433459) have to be considered fainter than R=20 (see Table 
2). Two sources (ELAISC15_J003034-433428 and ELAISC15_J003857-431719) were too close to bright stars 
and the optical identification was not possible. To estimate the identification reliability, we adopted the 
likelihood ratio technique as described by Sutherland & Saunders (1992) (see Ciliegi et al. (2003) for 
a detailed description). When possible counterparts fainter than R=23.0 were found, the identifications 
were not statistically discussed, and presented only as tentative. Taking into account all of the optical 
identifications, 332 ('--^82%) out of the 406 ELAIS-Sl sources have a reliable counterpart down to R'-^23 
(only 7 sources have 23.0<R<24.3), 2 are impossible to identify, and 72 (~18%) are to be considered fainter 
than R~23.0 (see Figure 3). In the following we will call these sources "EMPTY". 

3.3. K band and 1.4 GHz identifications 

The sample has been cross-correlated with the K-band catalog from the Two Micron All Sky Survey^ 
(2MASS), which reaches a K-band limit of ~15.6 (at 95% completeness level). In total 211 sources have been 
detected in the K-band. Since the whole SI field has been surveyed in the radio (1.4 GHz) with the Australia 
Telescope Compact Array to an average la level of 0.08 mJy (Gruppioni et al. 1999), at each ISOCAM 

position, we have searched for detection in the radio map down to 3 cr, finding 71 likely radio counterparts (61 
with measured redshifts; see Gruppioni et al. 2003 for details). The results of the multi-band identifications 
are reported in Table 2. 



4. The spectroscopic catalog 
4.1. Spectroscopic observations 

Spectroscopic observations of the optical counterparts of the ISOCAM SI sources were carried out at 

the 2dF/AAT and ESO Danish 1.5-m, 3.6-m and NTT telescopes during September and October 1998, 
September and November 2001, September and October 2002. A total of 15 ESO nights and 2 UK (from 
the Panel for the Allocation of Telescope Time; PATT) nights were allocated. 

The 2dF/AAT multifiber feed spectrograph was used with 200 fibers equipped with the 300B grism 
and the remaining 200 fibers with the 270R grism (about 10 A instrumental resolution). Because of poor 
weather conditions only half an hour of integration was possible during each night, reaching a magnitude 



®http:/ /www. nofs.navy.mil/. 

®See http;//www.ipac. caltech.edu/2mass. 
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limit of R~19.5. The ESO spectrographs (DFOSC, EF0SC2 and EMMI) were all used in single slit low 
resolution mode (10-20 A resolution) in the wavelength range 3500-9500 A, and reaching an average S/N 
ratio per resolution element of ~15-30 in the wavelength range 4000-7000 A. The reduction process used 
standard MIDAS and IRAF facilities. The raw data were sky-subtracted and corrected for pixel-to-pixel 
variations by division with a suitably normalized exposure of the spectrum of an incandescent source (flat 
field). Wavelength calibrations were carried out by comparison with exposures of He, Ar and Ne lamps. 
Relative flux calibration was carried out by observations of spectrophotometric standard stars (Oke 1990, 
Hamuy 1992, 1994). 

4.2. Classification 

We have spectroscopically classifled 290 (~87%) out of the 332 sources optically identifled (~71% of the 
whole 15 /um sample). The sample has been cross-correlated with the TYCH02 (H0g et al. 2000) and NED^° 
catalogs in order to identify known stars and extragalactic objects. Sixty-two TYCH02 stars were found. 
The remaining objects brighter than R=15 were visually inspected in the R band images and those suspected 
to be extended/extragalactic were spectroscopically observed. Fifteen sources with R<15 were considered 
to be bona fide stars on the basis of the visual inspection, while 14 were spectroscopically identified as stars 
at the telescope. In total we found 91 stars (all of them brighter than R=15) and 199 extragalactic sources. 

In order to classify the extragalactic sources we have first identified type 1 AGN (AGNl) when permitted 
lines with rest frame full width half maximum (FWHM) larger than 1200 Km/s were found. The remaining 
objects were classified according to two methods. 

The first method was used to identify type 2 AGN (AGN2), and to separate them from other emission 
line galaxies such as Hll/starburst and Low Ionization Narrow Emission Regions (LINER). All the objects 
with narrow 11/3 or [OIII]A5007 emission lines having EW (measured as the ratio between the integrated flux 
in the line, as resulting from a gaussian fit, and the flux of the continuum at the line wavelength) larger than 
'^3-5A (our line detection limit), were classified according to three diagnostic diagrams: 1) [OIIIJA5007/H/3 
ratio versus [NII]A6583/Ha ratio; 2) [OIII]A5007/H/3 ratio versus [SII]A6725/Hq! ratio; 3) [OIII]A5007/H/3 
ratio versus [OII]A3727/H/3 ratio ( e.g. VeUleux & Osterbrock 1987; Tresse et al. 1996; see Figure 4 for an 
example). The Ha and [NIIJA6583 lines were deblended using the IRAF task splot, which assumes Gaussian 
line profiles and a linear background. 

The above method is efficient in separating AGN from starburst activity, but does not provide infor- 
mation on the nature of a large fraction of our galaxies, which do not show [QUI] and H/3 emission lines. 
We have thus applied the classification scheme of Dressier et al. (1999) (see also Poggianti et al. 1999 
and Poggianti & Wu 2000) to all the "non AGN" galaxies. This method is based primarily on two lines, 
[OII]A3727 in emission and H(5 in absorption, which are good indicators of (respectively) current and recent 
star formation (SF) in distant galaxy spectra. Under this scheme the spectra have been divided in: 

1. e(a) (having 3A<EW([OII])<40A and EW(H(5)>4A), spectra of dust-enshrouded starburst galaxies; 

2. e(b) (having EW([OII])>40A), spectra with very strong emission lines of galaxies which are undergoing 
strong star formation; 



See http;//nedwww. ipac.caltech.edu/. 
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3. e(c) (having 3A<EW( [Oil] ) <40A and EW(H^) <4A) , typical spectra of spirals which have been forming 
stars in a continuous fashion; 

4. k (absent [Oil] and EW(H5)<3A), spectra of passive elliptical-like galaxies with neither ongoing nor 

recent star formation; 

5. k+a/a+k (absent [Oil] and EW(H(5)>3A), spectra of poststarburst galaxies with no current star for- 
mation which were forming stars at a vigorous rate in the recent past (last 1.5 Gyr). 

6. k(e), spectra similar to the k-type but with signs of at least one emission line. This class was introduced 
by us following Due et al. (2002). 

Following this classification, we found 25 AGNl, 23 AGN2, 37 c(a), 11 c(b), 68 e(c), 32 k(e) and 3 
k galaxies. No k+a/a+k galaxy was found. Thirteen objects could not be classified on the basis of this 
scheme because the wavelength range of [Oil] and B.5 lines was not observed (they were observed with the 
red grism of the 2dF spectrograph). Two of them were classified e(b) (starburst) galaxies according to the 
[OIII]A5007/H/3 and [NII]A6583/Ha ratios, while the remaining 11 galaxies were classified k(e) because no 
[OIII]A5007 and H/? line was detected, while the Ha emission line was present. In Figure 5 we show the 
composite spectra of each type of galaxy found under this classification scheme. 

We should stress that the number of identified AGN2 has to be considered a lower limit. XMM and 
Chandra surveys have revealed significant numbers of (mostly absorbed) AGN with L(2-10 keV)> 10^^ erg 
s""'^ (i.e. AGN2) whose optical spectra do not show sign of AGN activity (see e.g. Fiore et al. 2000; 2003). 
Recent (June- July 2003), 100 Ksec long, XMM observations of Sl-5 have revealed indeed that about 10-20% 
of the counterparts of our ISO sources classified as no- AGN galaxies on the basis of their optical spectra, 
can harbor an X-ray emitting AGN ( the BeppoSAX observations presented in Alexander et al. (2001) are 
too shallow to provide tighter AGN constraints). Under the reasonable assumption that all these objects 
are AGN2, the number of AGN2 in the ELAIS-Sl sample should be increased by at least a factor of 2. A 
more detailed discussion on the X-ray/MIR relationship of AGN and starburst galaxies is beyond the scope 
of this paper and will be presented by La Franca et al. (in preparation). 

In Table 2 we list the multiband associations to the 15 ^,m sources with the likelihood and reliability of 
the R band identifications, plus the redshift, the spectroscopic classification and the pLi^ luminosities at 15 
/xm and in the R-band. In Table 3 we present the line measurements of the spectra which include the calcium 
break (D4000; measured from the fluxes in the intervals 3750-3950 A and 4050-4250 A; see Dressier et al. 
1999) and the EW for the [OII]A3727, m, H/3, [OIII]A5007, Ha, [NII]A6583 and [SII]A6725 lines (positive 
values correspond to absorption lines). 

Figure 6 shows the distribution of the 211 ELAIS-Sl sources detected both in the R and K bands, in the 
R — K versus R — maglb plane, where we have defined magl5 = —2.5log{Si5{mJy)). This diagram shows 
a clear separation of the locus of the stars from the locus of the galaxies: for any given R — K color, stars 
populate a region with a two order of magnitude lower MIR to optical luminosity ratio. Note that only one 
early-type galaxy and two stars are located in the intermediate area of the diagram between those occupied 
by stars and galaxies. This means that stars can easily be separated from galaxies based purely on their 
infrared to optical flux ratio (see also Vaisanen et al. 2002). 
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4.3. Completeness 

The Sl_5 area is 99% spectroscopically complete (72/73) down to R=21.6. The Sl_rest area is 97% 
(216/222) complete down to R=20.5 (see Figures 7 and 8). These magnitude limits have been used for 
the statistical analysis of the evolution of AGNs and galaxies by Matute et al. (in preparation) and Pozzi 
et al. (2004) respectively. Both Figures 7 and 8 show a separation between stars and galaxies, the stars 
being optically brighter than galaxies, and also the appearance of a separate population of galaxies at faint 
magnitudes {R > 20). See §6.3 for a detailed discussion on this point. 

In the following we will call noIDs all the 116 ELAIS-Sl sources which are without spectroscopic 
identification (i.e., including also 72 EMPTY sources and 2 sources with no optical identification becouse 
too close to bright stars; see §3.2). 

5. Counts and contribution to the MIR background 

We have computed the contribution of each class of sources to the 15 /xm source counts in the mJy 
flux region. The integral counts for star-forming galaxies and AGNs (as defined on the basis of the optical 
spectra) are reported in table 4 (see Figure 9). As already discussed in section §4.2, the counts of AGN 
should be considered as lower limits. All ELAIS-Sl extragalactic sources (including the noIDs), with 15 /xm 
fluxes brighter than 0.6 mJy, produce a 15 iim background of i'I(i') = 0.34 nW m~^sr~'^. This corresponds 
to ^12% of the total extragalactic background estimated by Metcalfe et al. (2003) of z^I(z^)=2.7 nW m~^ 
sr^^. 

The sources spectroscopically identified with galaxies (no AGN) produce a flux of i/I(i/)=0.16 nW m~^ 
sr~^. The empty fields (mainly objects fainter than R~23) produce ul{u)=0.14: nW m~^ sr~^, while the 
AGNs produce z/I(z/)=0.04 nW m~^ sr"-'^. If we assume that the noID sources correspond mainly to starburst 
galaxies, we derive a '--^90% contribution of starburst galaxies and a ^^10% contribution of AGNs to the 15 
yum background at fluxes brighter than 0.6 mJy. If our preliminary estimate of the fraction of of hidden 
X-ray emitting AGN2 among the galaxies is taken into account, these figures may well be consistent with 
the higher estimates (~15-20%) for the contribution of AGN derived by Alexander et al. (2002), and Fadda 
et al. (2002). 

6. Nature of the 15 /^m extragalactic source population 

6.1. Redshifts and luminosities 

In Figure 10 the redshift distribution of the spectroscopically identified sources is shown. While AGNl 
reach redshifts up to 2~3, all AGN2 and galaxies have redshifts lower than 1, with an average of <2:>'~ 0.3. 
A possible hint of structure at 2; ~ 0.15 and 2; ~ 0.20 is present. 

We computed the vh^ luminosities of our sources at 15 /um (L15) and in the R-band (Lr)- For each 
class of objects a K-correction was applied according to the spectral classification (see Figure 11). In the 
infrared the K-corrections were based on: a) the SED from Elvis et al. (1994) for AGNl, b) the SED of 
NGC1068 for AGN2 (Sturm et al. 2000) and, c) the SED of M82 for all of the remaining sources, except the 
absorption line galaxies (for which M51 was considered). In the optical, the K-correction from Natali et al. 
(1997) was used for AGNl, while for AGN2 the K-corrections were computed by using the average of our 
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corresponding optical spectra. For the galaxies the SEDs of M82 and M51 were used, consistently with what 
has been done for the 15 /zm flux. The model SEDs of both M82 and M51 were taken from the GRASIL^^ 
library of models (Silva et al. 1998). For M82, the MIR (5 - 18 /um) has been replaced with the observed 
ISOCAM CVF spectrum (Forster Schreiber et al. 2003). 

In Figure 12 the L15 luminosities are shown as a function of redshift. A luminosity dependence of the 

classes of objects is evident. AGNl arc the most luminous and most distant sources with < logLi^>^ 12.1, 
while AGN2 are less luminous and less distant with < logLi^>^ 10.6. Finally, the least luminous and least 
distant sources are the galaxies with < logLi^>^ 9.9. 

6.2. The optical-MIR SED of starburst galaxies 

In Fi gure 13, Lir^/Ln versus £15 is shown for all of the ELAIS-Sl 15 /im galaxies plus sources from the 
spectroscopic identifications in the HDF-S from Franceschini et al. (2003), the HDF-N (ISOCAM data from 
Aussel et al. 1999; spectroscopic data from Cohen et al. 2000; magnitudes from Hogg et al. 2000), and the 
local sample of IRAS 12-/im galaxies of Rush, Malkan & Spinoglio (1993; RMS; after conversion to 15 /zm 
by using the M51 and M82 SED for normal and starburst galaxies respectively). A clear relation which is 
valid for all of the considered surveys is found, although it involves only objects with optical identification. 
The best fit to the data is 

log{Li5/Lii) = 0A7logLi5[LQ] - 5.02. 

The la dispersion of the relation is 0.32. Note that type 1 AGN have been excluded because of 
the different mechanism responsible for their luminosities, placing them on a different Lis/L/j versus L15 
plane location with respect to galaxies. On the contrary, type 2 AGN follow the same relation as all the 
other galaxies. This relation implies that the average optical-IR spectral energy distribution of galaxies is 
luminosity dependent, with more luminous galaxies having larger Li^/Lr ratios. A two order of magnitude 
increase in luminosity implies a 10 times larger Li^/Lr ratio. This relation is also discussed and taken into 
account for the analyses of the evolution of starburst galaxies and AGNs presented in two companion papers 
by Pozzi et al. (2004) and Matute et al. (in preparation). 

6.3. The nature of the sources fainter than the spectroscopic or photometric limits 

As already shown by Gruppioni et al. (2002), the 15-/im extragalactic counts show an upturn at fluxes 
fainter than ~1 my (see also Elbaz et al. 1999), which is interpreted as due to the contribution of the strong 
cosmological cvohition of the starburst population. This interpretation is discussed in detail by Pozzi et al. 
(2004). Anyway the counts (see Figure 9) show that the upturn is mostly due to the contribution of the 
noIDs sources (i.e. sources fainter than our spectroscopic or photometric limits), which represent about 
half of the extragalactic sources at fluxes ~0.6-l mJy. These sources could be either optically less luminous 
(absorbed?) galaxies at redshift similar to those of the already spectroscopically identified sample, or higher 
redshift sources. We favor the second hypothesis on the basis of the following discussion. 

We have tried to estimate a redshift for the noIDs sources in order to distinguish between the two 
different scenarios described above. We have combined our sample with the fainter HDF-S and HDF-N 



available at http://web.pd.astro.it/granato/grasil/modlib/modlib.html. 
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samples. Then we have tried to derive a flux-redshift relationship applicable to the ELAIS-Sl unidentified 
sources in order to estimate their redshifts. After several trials using multiple polynomial relations between 
log{z) and both 15 jim flux {logSi^) and magnitude {R), we have assigned a redshift to the 116 noIDS 
sources on the basis of their R magnitudes only {R = 23 was assumed for all the unidentified 76 {EMPTY) 
sources), by considering the following empirical relation and its spread: 

log{z) = -6.350 + 0.443 R - 0.007 + G(0, cr„.;)> 

where Urei (= 0.15) is the Icj dispersion of the relation and G(0, Urei) is a Gaussian distribution with center 
and width Orei (see Figure 14, Top). Note that the fit is non-linear in log{z) vs. R, since we found that 
a linear fit to the data extrapolated to faint magnitudes was providing too high redshifts, inconsistent with 
the k-correction limit of ISOCAM {z ^ 1.5) and with the median redshifts of i? = 24 — 25 objects from the 
HDF-S and HDF-N. Using the estimated redshifts, we have computed MIR and optical luminosities for the 
unidentified sources (by assuming all of them being star-forming galaxies, thus using M82 SED) and we have 
placed also these objects in the Li^/Lji vs. L15 diagram (see the bottom panel of Figure 14). 

According to our estimates, the unidentified objects in our survey are likely to be rather luminous MIR 

objects (_Li5 ranging between ~ 10^"'^ Lq and a few 10^^ L©, in the LIG/ULIG region), and their optical and 
MIR estimated luminosities appear to be consistent with (and extending to higher luminosities) the Li^/Ln 
vs. L15 correlation found for the identified sources. The estimated redshift distribution for the unidentified 
sources is shown in Figure 15, where it is compared with the distribution found for the spectroscopically 
identified sources and with the model fitting derived from our data (i.e. the 15 fim luminosity function in SI; 
Pozzi et al. 2004). The agreement between the model and the data (either observed and estimated) is good. 
It results that the unidentified sources are expected to fill the secondary high-z peak of the model-predicted 
redshift distribution. 

The double peaked redshift distribution (see Figure 15) is caused by the 15 /zm k-correction. As shown 
in Figure 11, at z=0.5 the 15 fiia observed fluxes decrease by ^30%. For our flux limited sample this implies 
that the ELAIS-Sl galaxies at 2;=0.5 correspond to 30% higher intrinsic luminosities. But at 2;=0.5 the 
ELAIS-Sl galaxies sample the bright steep slope of the luminosity function (log(Li5)>ll > L*), where 30% 
higher intrinsic luminosities correspond to an order of magnitude lower density of galaxies. This decrease of 
density of galaxies is also evident in the R-magnitude distribution (sec Figures 7 and 8), where, due to the 
k-correction effects, a lack of sources with R~20-21 (and z~0.5) is evident. The R-band histogram (Figure 
7) shows a triple peaked distribution: stars populate the brightest peak at R~10-12, 2;<0.5 galaxies populate 
the intermediate peak at R~17-19, and higher redshift (2;>0.5) galaxies populate the fainter peak that can 
be observed at R>21-22. 

Under the alternative hypothesis that ELAIS-Sl noIDs were instead lower redshift (^<0.5) sources 
suffering from dust extinction, it would be more appropriate to consider the 15-/L(m flux - z relation rather 
than the R — z relation (because affected by dust) to estimate their redshifts. In this case, we would obtain 
a redshift distribution for the noIDs sources similar to that of the optically identified, but they would be 
located in a separate region of the Lic,/Lr vs. _Li5 diagram (i.e. at similar _Li5, but with at least one 
order of magnitude higher Li^/Lji). This would imply the existence of a absolutely different population of 
objects, with much higher level of obscuration (without continuity) than the identified ones, but at similar 
flux density and redshifts. 

The first hypothesis is also supported by the faint spectroscopic identifications in the HDF-N (Aussel et 
al. 1999) and CFRS 1452+52 (Flores et al. 1999) fields which have mainly found a population of ^=0.4-1.0 
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galaxies. For all these reasons we conclude that the vast majority of the faint (R>20.5) ELAIS-Sl noIDs 
objects should belong to the same population as the identified sources, but at higher redshifts (2;=0.5-1.5), 
and that consequently the present sample is virtually spectroscopically complete up to z=Q.^, as far as the 
star forming galaxies are concerned. 

6.4. The nature of the MIR galaxies 

As we have seen, our sample is virtually complete up to z ~ 0.5. Up to this rcdshift we have identified 
167 extragalactic objects: 3 AGNl, 20 AGN2, 33 e(a), 9 e(b), 67 e(c), 32 k(e), and 3 k galaxies. No k+a/a+k 
galaxy was found. The average redshift of the sample is <z>=Q.2Q. The number of e(a) galaxies rises to 
37 if we classify also the AGN2 according to the classification scheme of Dressier et al. (1999) (see §4.2). 
The majority of the sources with 2;<0.5 is composed by normal star forming spirals (e(c), 41%), and by 
absorbed star forming galaxies (e(a), 20%; 23% including the AGN2). The total fraction of optical spectra 
which show emission lines which are indicators of star formation (e(b)+e(c)+e(a)+k(e)) is 84%. 14% of the 
sources arc AGNs, and only 2% of the sources do not show emission lines (k galaxies). This confirms that 
the MIR and the optical data agree on revealing both ongoing star formation and accretion driven (AGN) 
activities. The lack of k+a/a+k galaxies in the ELAIS-Sl sample suggest that these galaxies are not strong 
MIR emitters. This result is consistent with the lack of emission lines in their optical spectrum and with 
their post-starburst interpretation (see Poggianti and Wu 2002). 

The c(a) galaxies are believed to be associated with dust enshrouded starburst galaxies. Following 
Poggianti et al. (1999) and Poggianti & Wu (2000), the spectra of e(a) galaxies are explained as the result 
of selective dust extinction which affects the youngest, most massive, stars in HII regions much more than 
the older stellar population responsible for the continuum flux. A search for similar spectra in the local 
universe (Poggianti et al. 1999) revealed that they are frequent among merging/strongly interacting galaxies 
(~40%), while they are scarce in optically selected surveys of field galaxies at low z: 7-8% in the Las 
Campanas Redshift Survey and in Kennicut (1992). Dressier et al. (1999) and Poggianti et al. (1999) found 
that e(a) galaxies constitute about 10% of their sample of cluster and field galaxies at 2:~0.4-0.5. 

At 15 /um most (71%) of the sources with optical spectroscopy in the CFRS are classified as e(a) galaxies 
and have a mean redshift <z>=0.76, while at 2;=0.18, in the cluster of galaxies Abell 1869, Due et al. (2002) 
find 25% of 15 /im galaxies classified as c(a). 

The observed 20-23% fraction of e(a) galaxies in the ELAIS-Sl sample at <2;>=0.20, is thus higher 
than observed in optically selected samples and in agreement with the estimates at similar redshifts by Due 
et al. (2002) in a MIR selected sample. This confirms that e(a) galaxies are dusty absorbed star forming 
galaxies which are preferentially selected at longer wavelengths (MIR), and whose frequency increases with 
increasing redshift, as the cosmic star formation rate. 

We have compared the R — K colors of our sources to those of the field galaxies and we did not find 
evidences of reddening. Becouse of the bright K-band limit magnitude we limited our comparison in the 
magnitude range R=15-17.5, and consequently the statistic are poor. The average R-K color is 2.04 and 
2.27 for FLATS and field galaxies respectively, with a Icr dispersion of ^0.50 for both samples. On the 
contrary, dust extinction seems to affect the spectral lines of our objects. In figure 16 the EW of the [Oil] 
line is compared to the EW of the Ha+[NII] line. The EW{[OII])/EW{Ha + [Nil]) ratio is related to the 
differential extinction due to dust, the [Oil] emission being more affected than Ha emission due to its shorter 
wavelength. The dashed line in Figure 16 is the relation found by Kennicutt (1992) for local field galaxies 
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{EW{[OII])/ EW{Ha + [NII])=OA). Most of our sources lie below that relation, and the presence of upper 
limits (i.e. [Oil] non-detections) would even lower the average value of the EW{[OII])/EW{Ha + [Nil]), 
implying a larger amount of extinction in our data than in the local sample of Kennicutt (1992). The average 
ratio for the galaxies having both [Oil] and Ha+[NII] lines is 0.36±0.15, 0.35±0.03, 0.29±0.02 for e(b), e(c) 
and e(a) galaxies respectively. Thus, as expected, at least for the normal star forming (e(c)) and absorbed 
star forming (e(a)) galaxies the difference from the local observed ratio (0.4) is significant. 

In order to disentangle a possible dependence on redshift or luminosity of the average extinction of the 
galaxies, we have combined our sample with a subsample of 39 galaxies from the local sample of RMS for 
which EW measures for [Oil] and Ha-|-[NII] were available in the literature (see Alexander & Aussel (2000) 
for a discussion of this sample) . The data show that the extinction affecting the lines is an increasing function 
of luminosity while a correlation with redshift is not statistically significant: a partial correlation Kendall 
test gives in fact a 10~^ probability that the EW is not dependent on luminosity, and a 0.04 probability 
that the EW is not dependent on redshift. Figure 17 (top) shows the distribution of the RMS and ELAIS 
galaxies in the luminosity-rcdshift plane. We have then selected the galaxies with 0.010<z<0.055 in order 
to show the dependence on luminosity (Figure 17, middle), and the galaxies in a small luminosity interval 
(9</o5'Z/i5(i©)<10) in order to show the not significant trend with the redshift (Figure 17, bottom). These 
results suggest an increase of dust extinction with luminosity. Similar results, based on a smaller sample, 
were found also by Pozzi et al. (2003) in one of the smaller but deeper ELAIS fields, S2. 

A quantitative estimate of the amount of optical extinction affecting the ELAIS-Sl MIR galaxy sample 
can be obtained by the measure of the observed Balmer decrement in the Ha and H/3 lines. We have made 
this measure on two composite spectra for each class of galaxies: one including only spectra with detected 
Ha and H/3 emission lines, and one including also spectra with no detected H/3 line. Spectra showing H/3 
lines in absorption were excluded. This method was adopted in order to minimize the effects of underling 
H/3 absorption which causes the observed Ha/H/3 ratio to be larger. For this reason the extinction estimates 
we derive have to be considered upper limits. In Table 5 our estimates in the case of B recombination with 
a density of 100 cm~^ and a temperature of 10000°K (Osterbrock 1989) are shown. The higher values of 
E(B-V) were always obtained when the spectra without H/? detection were included in the composite spectra. 
These values correspond to an extinction Ay=3.2x£'(fi — V) of about 1.5, 3.5, 4.0, 4.5 mag for e(b), e(c), 
e(a) and k(e) galaxies respectively. These value refers to the absorption of the lines emission; typically the 
absorption on the contimmm is a factor of two smaller (see also Pozzi ct al. 2003, and Poggianti & Wu 2000 
for similar results). As expected, the larger extinction values are found for e(a) and k(e) galaxies, reinforcing 
the interpretation that the e(a) signature is associated with dust extincted starbursts. 



6.5. The star formation rate 

We have estimated the star formation rate (SFR) of the ELAIS-Sl sources from their 15 /um infrared 
luminosity. The SFR estimator has been computed as: 

SFR[MQ/yr] = 11.25 x lQ-^^Li5[WHz-\ 

assuming a Salpeter initial mass function (IMF, Salpeter 1955) between 0.1 and 100 Mq, a dependence 
of SFR from Leo;:^™ as given by Mann et al. (2002; SFR[MQ/yr] = 0.9 x lO-'^^Leo^,m[WHz-'^]), and a mean 
value of the ratio < Li5f^rn[WHz^^]/LQ0fim[WHz~^] >~ 0.08, typical of starburst galaxies such as M82. 
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In Figure 18 the SFRs as a function of redshift are shown. Type 1 AGN have been excluded as their MIR 
luminosities are not attributed to star forming activity, but to reprocessed higher frequency luminosity driven 
by the accretion on a super-massive black hole (10^-10^° M©). For the sake of comparison the estimates 
from the spectroscopic identifications in the HDF-S, and the HDF-N are shown. The ISOCAM observations 
of the HDF-S and HDF-N reach 15 /im fluxes down to ^^0.1 mJy (i.e., more than six times fainter than 
our survey) and probe a different region of the SFR vs. z plane. The shallower wide-area ELAIS survey is 
complementary to the deep surveys and is required to make a more accurate estimate of the SFR history 
of the Universe (see Pozzi et al. 2004). The ELAIS-Sl sources have mean SFRs of 13 Mg/yr, 45 Mg/yr 
and 242 Mg/yr, in the redshift ranges z<0.2. 0.2<z<0.4 and z>0.4 respectively. Obviously this should not 
be attributed to an evolutionary trend, but to the luminosity-redshift relation which is present in any flux 
hmited sample like ELAIS-Sl. 

We have found a good correlation between the SFRs derived from the 15 /xm luminosity and those 
derived from the 1.4-GHz luminosity (computed as SFB\^^\ = 1.1 x W~^^{qjj^)°'Li,4ghz[-^], assuming a 
Salpctcr IMF between 0.1 and 100 Mq, and a power-law spectrum with a spectral index a ^0.7; Su(xi^~°'). 
The correlation between the SFRs obtained through the two different indicators reflects the radio-MIR 
luminosity correlation found for ELAIS-Sl galaxies by Gruppioni et al. (2003; see also Elbaz et al. 2002). 

7. Summary 

We have presented multiband (R, K, 1.4 GHz) and optical spectroscopic identifications of a reliable 
subsample of 406 15-/Lim sources from the ELAIS-Sl survey over the flux density range 0.5< Si^fj^m <150 
mJy. These observations substantially contribute to the understanding of the evolution of star-forming 
galaxies and AGN in the MIR, which is presented in two companion papers by Pozzi et al. (2004) and 
Matute et al. (in preparation) respectively. This paper has described the details of the ground-based follow- 
up observations and data reduction, and has presented optical identifications for ~82% of the 15 fim sources 
down to R'~23, and spectroscopic identification for ~88% of the optically identified sample. The areal 
coverage and completeness of the survey has been described in detail, in order to allow statistical analyses. 
The source counts for star-forming galaxies and AGN have been derived: the ELAIS extragalactic sources 
down to Si5^m=0.6 mJy have a density of 370±14/deg^, and produce a 15 /xm background of j^I(j^) = 0.34 
nW m~^sr~^. This corresponds to ~12% of the total extragalactic background estimated by Metcalfe et al. 
(2003). Taking into account the optical based classification, we estimate an AGN contribution of at least 
-^6% to the MIR sources brighter than Si5^m=0.6 mJy. A larger fraction is expected if the classification 
will be based on the X-ray luminosity. The line ratios of the optical spectra of the MIR galaxies suggest 
that they are affected by substantial reddening, whith the more MIR-luminous sources having larger dust 
extinction. A general trend is foimd in the MIR-optical SEDs of IRAS and ISO galaxies whith the more 
MIR-luminous sources having larger Li^/Lji ratios. According to our analyses we conclude that most of 
the ELAIS galaxies fainter than R~23 should have redshift greater than 0.5 and are likely to be LIG-ULIG 
galaxies. As a consequence, the present sample of star-forming galaxies has to be considered virtually 
complete up to .2 = 0.5. 

Based on observations collected at the European Southern Observatory, Chile, ESO N°: 57.A-0752, 58. B- 
0511, 59.B-0423, 61.B-0146, 62.P-0457, 67.A-0092(A), 68.A-0259(A), 69.A-0538(A), 70.A-0362(A); PATT 
proposal P076). This paper is based on observations with the Infrared Space Observatory (ISO). ISO is 
ane ESA project with instruments funded by ESA Member States (especially the PI countries: France, 
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Germany, the Netherlands and the United Kingdom) and with the partecipation of ISA and NASA. This 
research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics 
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Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and 
the National Science Foundation. This research has been partially supported by ASI contracts ARS-99-75, 
ASI I/R/107/00, ASI I/R/113/01, MURST grants Cofin-98-02-32, Cofin-00-02-36. IM acknowledges a PhD 
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Fig. 1.— ELAIS-Sl field divided into the 9 ISOCAM rasters. The total area is 4.2 deg^ The numbers refer 
to the individual rasters. The 168 CCD pointings and the 406 ELAIS-Sl sources are shown by dotted lines 
and filled points, respectively. 
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Fig. 2. — Completeness as a function of the 15 /um flux of the Sl-5 (continuous line) and the Sl-rest (dashed 
line) regions. 
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Fig. 3. — 15/im flux histograms of the 406 ELAIS-Sl sources in the Sl-5 [Top) and Sl-rest [Bottom) 
regions. The shaded histograms are the opticaUy identified sources, and the cross-shaded histograms are the 
spectroscopically identified sources. 
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Fig. 4. — Diagnostic diagram for the classification of emission line galaxies. The dashed lines separate the 
regions of AGN2, LINERs and starburst galaxies (see e.g. Veilleux & Osterbrock 1987 and Tresse et al. 
1996). The errorbars correspond to 10% uncertainties. 
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Fig. 5. — Composite spectra of the various type of galaxies classified in the ELAIS-Sl survey (see text). 
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Fig. 6. — R — K color versus R — maglh color diagram of the 211 ELAIS-Sl sources detected both in the R 
and K bands, maglb is a 15 /xm magnitude computed as magl5 = —2.blog{S\^{mJy)) . 
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Fig. 7. — R-band magnitude histograms of the optically identified 15 fim sources in the Sl-rest {Top) and Sl- 
5 (Bottom) regions. The shaded histograms are the R-band distributions of the spectroscopically identified 
sources. 
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Fig. 8. — ELAIS-Sl spectroscopically identified sources in the Si5^m-R-band plane. Upper panel: sources 
in the Sl-5 area. Lower panel: sources in the Sl-rest area. The dashed hnes are the optical limits of the 
spectroscopically complete samples (see §4.3), while the upper arrows are sources without optical counterpart 
on CCDs. 
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Fig. 9. — Integral counts from the ELAIS-Sl survey for a) all the extragalactic sources, b) the galaxies plus 
the no identified sources which are assumed to be galaxies, c) the AGNs and d) the no identified sources. 
The counts of galaxies and AGN corrected for the contribution of X-ray emitting AGNs are also shown. The 
no identified sources {noIDs) includes 72+2 EMPTY sources plus 40 (mainly faint: 20.5<R<23.3) sources 
without redshift. Errors are at the la confidence level, based on Poisson statistic. 
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Fig. 10. — Redshift distributions of the spectroscopically identified sources in ELAIS-Sl. Top: AGNI; 
Middle: AGN2; Bottom: all the extragalactic sources which are not AGNI. 




Fig. 11.— K-correction ( log(S(z)/S(0)) ) in the IS0CAM-LW3 15 /xm and R-band filters for type 1 AGN, 
type 2 AGN and Starburst galaxies. For Starburst galaxies the MIR K-correction with M82 is shown as it 
is used in the vast majority of the cases (see §6.1). 



- 28 - 



0.1 




o o 



I ' ' 

O CP 

O O qO 



— # ^ o 

^ • • 



+ 



A 
^ A 



+ A 



_l I I I I I I I L 



O AGNl 

• AGN2 

A e(a) 

□ e(b) 

+ e(c) 

A k(e) 

V k 



J I I I I I I I I I L 



8 9 10 11 12 13 

Fig. 12. — Luminosity at 15 /im of the ELAIS-Sl spectroscopically identified sources as a function of redshift. 
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Fig. 13. — Lis/iii ratio as a function of L15 for different MIR surveys. Filled circles are the ELAIS-Sl 

sample; open squares and asterisks are data from HDF-S and HDF-N, respectively; filled triangles are IRAS 
galaxies from Rush, Malkan and Spinoglio (1993). Type 1 AGN have been excluded, due to the different 
mechanism responsible for their luminosities. 
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Fig. 14. Top: R — z plane with the lower limits estimates of the redshift for the noIDs sources as 
described in §6.3. Filled circles are ELAIS-Sl data points; filled squares, asterisks and filled triangles are 
data from HDF-S, HDF-N and RMS respectively. ELAIS-Sl sources with optical counterparts, but lacking 
redshift information are shown as open circles, while empty fields are represented as pointed- up arrows (all 
at i? = 23). The dashed line is the best-fitting polynomial function used for estimating redshifts for sources 
without spectroscopic identification. Bottom: Same as in figure 13, but with also the lower limits for the 
noIDs sources shown. Symbols are the same as in the top panel. Data from RMS are not shown in this plot 
for greater clarity. 
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Fig. 15. — Measured redshift distribution for the spectroscopically identified sources in SI (solid histogram) 

and estimated redshift distribution for the unidentified ones (dashed histogram). Data are compared with 
the prediction of the model derived from the 15 /xm luminosity function obtained from our data and fitting 
the source counts (see Pozzi et al. 2004), shown as dotted line. 
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Fig. 16.— EW of [Oil] versus EW of Hq!+[NII] for ELAIS-Sl sources (symbols as in Figure 6), and 39 IRAS 
galaxies from Rush, Malkan & Spinoglio (1993, small dots), using data from literature;. The pointed-down 
arrows are sources with Ha+[NII] measurement but no [Oil] detected. The dashed lino shows the relation 
EW{[OII]) = OAEW{Ha + [Nil]) found by Kennicutt (1992) for local field galaxies (see §6.4). 
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Fig. 17. — Top - Luminosity-rcdshift distribution of the EL AIS-Sl galaxies (filled circles) plus 39 lists' galax- 
ies from RMS (triangles). Dashed lines show the regions selected for disentangling the luminosity /redshift 
dependence of the reddening (see §6.4). Middle - EW{[OII])/EW{Ha+[NH]) ratio versus 15 fim luminos- 
ity for the galaxies with 0.010<2;<0.055 (symbols as above). Bottom - EW{[OII])/EW{Ha + [Nil]) ratio 
versus redshift for the galaxies with 9 < logLi;,{LQ) < 10 (symbols as above). The remaining ELAIS-Sl 
plus RMS galaxies are represented with small filled circles. 
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Fig. 18. — Rates of tlie star formation based on the estimates of tiie MIR (15 /um) luminosity, as a function 
of redshift. Type 1 AGN have been excluded as their MIR luminosity is not attributed to the star forming 
activity. 



Table 1. The completeness function of Sl-5 and Sl-rest 



Flux Sl-5 Sl-rcst 






.45 


0. 


.001 


0.000 





.50 


0, 


.007 


0.000 





.55 





.023 


0.000 





.60 


0, 


.053 


0.000 





.65 





.110 


0.001 


0, 


.70 


0, 


.167 


0.004 





.75 


0. 


.239 


0.009 





.80 


0. 


.312 


0.019 





.85 


0. 


.379 


0.031 





.90 


0. 


,442 


0.049 





.95 


0. 


.501 


0.076 


1. 


.00 


0. 


.545 


0.107 


1 


.05 


0, 


.585 


0.147 


1, 


.10 


0, 


.625 


0.189 


1, 


.15 


0, 


.668 


0.232 


1 


.20 


0, 


.708 


0.278 


1 


.25 


0, 


.743 


0.323 


1 


.30 


0, 


.776 


0.373 


1 


.35 


0. 


.807 


0.415 


1 


.40 


0. 


.841 


0.465 


1 


.45 


0. 


.866 


0.506 


1 


.50 


0. 


.889 


0.559 


1 


.60 


0. 


,929 


0.646 


1 


.70 


0, 


.957 


0.699 


1 


.80 


0, 


.977 


0.743 


1 


.90 


0, 


.991 


0.781 


2 


.00 


1, 


.000 


0.813 


2 


.10 


1, 


.000 


0.841 


2 


.20 


1. 


,000 


0.867 


2 


.30 


1. 


.000 


0.888 


2 


.40 


1. 


.000 


0.906 


2 


.50 


1. 


,000 


0.921 


2 


.60 


1. 


.000 


0.934 


2 


.70 


1. 


.000 


0.945 


2 


.80 


1, 


.000 


0.954 


2 


.90 


1, 


.000 


0.962 


3 


.00 


1, 


.000 


0.969 


3 


.10 


1. 


.000 


0.974 


3 


.20 


1. 


.000 


0.979 


3 


.30 


1. 


.000 


0.983 


3 


.40 


1. 


.000 


0.986 


3 


.50 


1. 


.000 


0.989 


3 


.60 


1, 


.000 


0.991 


3 


.70 


1, 


.000 


0.993 


3 


.80 


1, 


.000 


0.995 


3 


.90 


1, 


.000 


0.996 


4 


.00 


1, 


.000 


0.997 


4 


.10 


1, 


.000 


0.999 


4 


.20 


1, 


.000 


1.000 



Table 2. Sources in the ELAIS SI Sample 



Name 


# 


Fl5 


ao(J2000) 


Sq (J2000) 


Ott 


R 


K 


Lik 


Rcl 


Fl.4 




CI 


Ll5 




Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(S) 


(!)) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


ELAISC15_J003957-432013 


S 


3.01 


00:39:57.85 


-43:20:13.4 


1.2 


17.21 


15.32 


363.16 


0.999 


< 0.24 


0.128 


5 


9.78 


10.27 




ELAISC15_J003958-441511 


9 


25.39 


00:39:58.67 


-44:15:10.6 


2.0 


10.41 


7.39 


999.99 


1.000 


< 0.30 


0.000 


8 








ELAISC15_J004000-431325 


8 


1.01 





























ELAISC15_J004004-433710 


9 


1.24 





























ELAISC15_J004007-432516 


8 


3.11 


00:40:06.88 


-43:25:18.2 


2.5 


12.21 


9.46 


999.99 


1.000 


< 0.30 


0.000 


8 








ELAISC15_J004009-434424 


9 


2.62 


00:40:09.26 


-43:44:25.3 


0.9 


19.13 




97.33 


0.997 


< 0.30 


0.188 


2 


10.15 


9.87 




ELAISC15_J004010-432012 


8 


1.03 






























Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. The columns are as follows: (1) Name of 
the 15-/jm source as in Lari et al. (2001). (2) Raster number. (3) The 15-/Lim flux in mjy units. (4-5) The ra and dec of the optical counterpart, J2000. (6) Offset of the optical 

counterpart from the 15-fim source position in arcsec units. (7-8) The R and K band Vega magnitudes. (9) Likelihood of the optical identification; 999.99 was assigned to stars (see 
§3.2). (10) Reliability of the optical identification (see §3.2). (11) 1.4 GHz flux or 3cr upper limit in m.Ty units. (12) Redshift. (13) Spectroscopic class where: means no spectra 
available, 1 is AGNl, 2 is AGN2, 3 is e(a) galaxy, 4 is e(b) galaxy, 5 is e(c) galaxy, 6 is k(c) galaxy, 7 is k galaxy and S is star. An L is appended for objects classified LINERS 
according to the line ratios diagnostics (see §4.2). (14) log(iyii^) IS-^im luminosity in solar units. (15) logivLf) R band luminosity in solar units. (16) Notes, O is for outside CCD 
observations, M is for multiple optical counterparts, and B for Broad Absorption Line (BAL) QSO. 



Table 3. Rest frame line measurements 



Name 


Fl5 




CI 


D4OOO 


EW([OII]) 


EW(H5) 


EW(H/3) 


EW([OIII]) 


EW(Ha) 


EW([NII]) 


EW([SII]) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



ELAISC15_.T002904-432415 


1.61 


0.207 


5 


1.20 


-8.9 


3.1 


-3.7 


0.0 


-56.7 


-10.6 


-12.0 


ELAISC15_J0029 15-430333 


2.04 


0.417 


5 


1.23 


-39.2 


3.4 


-10.3 


-13.9 








ELAISC15_J002924-432233 


2.34 


0.374 


5L 


1.17 


-39.4 


0.0 


-3.4 


0.0 


-36.4 


-47.3 


-46.3 


BLAISC15_J002939-43n625 


3.73 


0.071 


5 


1.22 


-11.8 


3.6 


-1.9 


0.0 


-22.7 


-8.6 


-8.4 


ELAISC15.J002949-430703 


1.33 


0.302 


2 


1.40 


-7.4 


0.0 


0.0 


-25.0 


-15.8 


-23.8 


-12.9 


EL AISC15_J00300 1-432202 


1.47 


0.274 


5 


1.33 


-3.4 


3.4 


0.0 


0.0 


-28.1 


-8.0 


0.0 


ELAISC15_J00301 1-432947 


1.35 


0.084 


6 


1.S7 


-2.6 


0.0 


-1.1 


0.0 


-4.8 


-3.4 


0.0 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. The columns are as follows: 
(1) Name of the 15-/:tm source as in Lari et al. (2001). (2) The 15-/jrm flux in mjy units. (3) Redshift. (4) Spectroscopic class as in Table 2. (5) Calcium break. 
(6-12) Rest frame EW in A units of the [OII]A3727, H5, H/3, [Oin]A5007, Ha, [NIIJA6583, and [SIIJA6725 lines; positive values mean absorption. 
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Table 4. 


The integral 


extragalactic counts at 15 /xm 


Sl5 


All 


Galaxies 


AGN1 + AGN2^ 


AGNl 


AGN2^ 


(mJy) 


(N(>S)/dog2) 


(N(>S)/dog2) 


(N(>S)/dcE2) 


(N(>S)/dog2) 


(N(>S)/dcg2) 


0.60 


370±14 


349 ±13 


21.5±3.3 


13.6 ±2.6 


7.9 ±2.0 


0.76 


242±11 


221 ±10 


21.5±3.3 


13.6 ±2.6 


7.9 ±2.0 


0.95 


140.8±8.4 


122.2 ±7.8 


18.6±3.1 


10.7 ±2.3 


7.9 ±2.0 


1.20 


70.1±5.9 


57.6 ±5.4 


12.5±2.5 


6.1 ±1.7 


6.4 ±1.8 


1.51 


42.4±4.6 


32.7 ±4.0 


9.7±2.2 


4.8 ±1.6 


4.9 ±1.6 


1.90 


22.8±3.7 


20.7 ±3.2 


7.1±1.9 


3.2 ±1.3 


3.9 ±1.4 


2.39 


19.3±3.1 


13.3 ±2.6 


6.0±1.7 


2.9 ±1.2 


3.1 ±1.2 


3.01 


14.1±2.7 


9.3 ±2.2 


4.8±1.5 


1.9 ±1.0 


2.9 ±1.2 


3.79 


9.5±2.2 


5.9 ±1.7 


3.6±1.3 


1.2 ±0.8 


2.4 ±1.1 


4.77 


7.3±1.9 


4.3 ±1.5 


3.1±1.2 


1.0 ±0.7 


2.1 ±1.0 


6.00 


5.7±1.7 


2.8 ±1.2 


2.8±1.2 


1.0 ±0.7 


1.8 ±1.0 


7.55 


3.6±1.3 


1.90±0.97 


1.7±0.9 




1.4 ±0.8 


9.51 


2.6±1.1 


1.19±0.77 


1.4±0.8 




1.2 ±0.8 


11.97 


1.9±1.0 


0.95±0.69 


1.0±0.7 






15.07 


1.19±0.77 


0.71±0.60 








18,98 


0.71±0.60 


0.48±0.49 









^Including sources without spectroscopic classification (nol Ds) . 
'^These values should be considered lower limits (see §4.2 and §5). 



Table 5. Mean E(B-V) 



Class E(B-V) 



e(b).... 


0.43- 


0.52 


e(c) 


0.82- 


1.59 


e(a).... 


1.02- 


1.42 


k(e).... 


0.64- 


1.87 



